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ABSTRACT 

Context. We have started a campaign to identify massive star clusters inside bright molecular bubbles towards the Galactic Center. The 
CN15/16/17 molecular complex is the first example of our study. The region is characterized by the presence of two young clusters, 
DB 10 and DB 1 1 , visible in the near infrared, an ultra-compact H n region identified in the radio, several young stellar objects visible 
in the mid infrared, a bright diffuse nebulosity at 8ii.m coming from PAHs and sub-mm continuum emission revealing the presence of 
cold dust. 

Aims. Given its position on the sky (Z = 0:58, b = -0?85) and its kinematic distance of 7.5 kpc, the region was thought to be a very 
massive site of star formation in proximity of the Central Molecular Zone. One of the two identified clusters, DB 1 1 , was estimated 
to be as massive as ~ 10 4 Mq. However the region's properties were known only through photometry and its kinematic distance was 
very uncertain given its location at the tangential point. We aimed at better characterizing the region and assess whether it could be a 
site of massive star formation located close to the Galactic Center. 

Methods. We have obtained NTT/SofI deep JHKs photometry and long slit K band spectroscopy of the brightest members. We have 
additionally collected data in the radio, sub-mm and mid infrared, resulting in a quite different picture of the region. 
Results. We have confirmed the presence of massive early B type stars and have derived a spectro-photometric distance of ~ 1.2 
kpc, much smaller than the estimated kinematic distance. Adopting this distance we obtain clusters masses of M DB io ~ 170Mq and 
Mdbii ~ 275Mq. This is consistent with the absence of any O star, confirmed by the excitation/ionization status of the nebula. No 
He i diffuse emission is detected in our spectroscopic observations at 2.113irm, which would be expected if the region was hosting 
more massive stars. Radio continuum measurements are also consistent with the region hosting at most early B stars. 

Key words. ISM: individual: CN15/16/17 - ISM: bubbles - Open clusters and associations: individual: DB 10, DB 1 1 - H n regions 

1. Introduction These differences highlight the need for further observational 

constraints on the star formation scenario close to the GC. 
Due to the high degree of interstellar extinction, very few young 

clusters are known in the Galactic Plane at distances larger than 2 

^ kpc from the Sun. Nevertheless, the central region of the Galaxy, Some studies suggest that several hundred clusters are ex- 

with its unique physical conditions, represent s a very interesting pected to reside along the line of sight towards the inner Galaxy 

laboratory for testing star formation theories dMorris & Serabvnl (Portegies Zwart et al. 2001 ), but very few of these clusters have 

1996). Recent studies give somewhat controversial results on the so far been identified. Unfortunately, if we look towards the cen- 

star formation efficiency (SFE) in the proximity of the Galactic ter of the Galaxy, the identification of clusters is complicated 

Cente r (GC): some suggest a SFE similar to the Milky Way by the severe crowding of fore/background stars. Spurious clus- 

Disk (lYusef-Zadeh et alJl2009h . others imply a reduced SFE to- ter detections are also possible due to the patchy nature of the 

wards the GC, based on a Galactic-wide comparison of dense gas interstellar medium, which may cause strong spatially varying 

tracers and active star formation tracers uLongmore et al.lE012l) . foreground extinction and consequently a variation in the star 

counts. Due to the strong extinction along lines of sight within 

* Based on observations performed at ESO's La Silla-Paranal obser- the Galactic Plane, the detection of distant clusters at optical 

vatory. Programme ID 085.D-0780. wavelengths is challenging, even for the most massive clusters 

** Member of the "International Max Planck Research School for located towards the center of the Milky Way. The best way to 

Astronomy and Cosmic Physics at the University of Heidelberg" identify new clusters in the inner Galaxy is therefore to use in- 

(IMPRS-HD). frared wavelengths, which are less affected by extinction. 
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Fig. 1. Left: Sofl J (blue), H (green) and K s (red) composite image of the CN15/16/17 complex. Right: Spitzer/IRAC 3.6^m (blue), 
5.8 p.m (green) and 8.0ixm (red) composite image for the same region. The white contours are NVSS radio continuum maps. The 
physical length-scale of 0.2pc is estimated at a distance of 1200 pc (see Sect.|5j- 



Over the past decade, the search for new clusters has 
gained fresh interest thanks to n ear infrared surveys such 
as D ENIS dEpchtein et all 1 19991) 2MA SS dSkrutskie et alJ 
120061). UKIDSS-GPS llLucas et al] 12008) and VISTA-VVV 
dArnaboldi et al.ll2007l: lMinnitietal]|2010l) . In parallel to these 
near infrared ground based surveys, mid infrared surveys of the 
Galactic Plane have been carried out using the IRAC camera 
on bo ard the Spitzer Space Telescope. IChurchwell et al.l ([2006, 
120071) found almost 600 molecular bubbles in the GLIMPSE I 
and II surveys. Molecular bubbles associated with H n regions 
are a tracer of young, massive clusters hosting O-type or early 
B-type stars. Within the sample identified in the GLIMPSE II 
survey (|/| < 10°, \b\ < 1°), 29 bubbles are associated with H n 
regions, meaning that they very likely host massive stars. 

We recently started a campaign aimed at characterizing the 
stellar content of these 29 regions, using a combination of imag- 
ing and long slit spectroscopy. In this paper we present the re- 
sults for CN15/16/17. We observed the region using the Sofl 
near infrared instrument mounted on the ESO-NTT telescope in 
La Silla, Chile. The observing strategy consisted of a combina- 
tion of deep JHKs imaging and long slit K band spectroscopy 
of the brightest candidate members. By combining photometry 
and spectroscopy we have been able to confirm the presence of 
early B type stars and to constrain the region's distance using 
their spectral type classification (see Sect. [5}- 

The CN15/16/17 complex of molecular bubbles is a star 
forming region (SFR) hosting young stars in different evolu- 
tionary phases (see Fig. [TJ. The region is projected towards 
the Galac tic Center (I = 0°58, b = -0°85) and was first de- 
tected by IChurchwell et al.l d2007l) . by visually searching the 
inner 20° of the Gala xy, using mid infrared data from the 
GLIM PSE II survey dBeniamin et all 120031: IChurchwell et alJ 
2009). The Spitzer/IRAC images of the region show a very 
pronounced diffuse emission in the 8ii.m channel, originating 
from PAH emission. Two stellar clusters are associated with 
the region and are visible in the near infrared. The clusters 



were first identified by iDutra & Bical (120001) using 2MASS im- 
ages (Skrutskie et al. 2006). One of them (DB11) has already 
emerged from its parental cloud and therefore its stellar popula- 
tion is detectable in the near infrared. The second (DB10) is still 
deeply embedded and the high extinction only allows the detec- 
tion of the brightest sources. Dutra et al. (2003) further studied 
the clusters using H and K$ band imaging also obtained with 
NTT/Sofl. In addition, a third, very deeply embedded SFR is 
present. This youngest region consists of a group of young stel- 
lar objects (YSOs) visible in the Spitzer images and corresponds 
to the IRAS 17470-2853 source. It is associated with a radio 
detecte d Ultra Compact H n region as well as several methanol 
masers dWalsh et al.lll998l) . 

In Fig.Q] left panel, we show a JHK$ composite image of the 
region from our Sofl observations. Objects DB10 and DB1 1 are 
clearly visible. DB11 is the central, larger cluster, while DB10 
is the smaller cluster west of DB 1 1 . The Spitzer/IRAC compos- 
ite image, right panel, is dominated by bright PAHs emission 
in the 8|xm channel. The contours indicate the radio continuum 
flux at 1.4 GHz fro m the NRAO VLA Sky Survey (NVSS, see 
ICondon et al.| [l998f). From the radio contours it is possible to 
identify two H n regions, one associated with DB11, and the 
other with the aforementioned Ultra Compact H n region, corre- 
sponding to the third, deeply embedded, SFR in the complex. In 
the latter, a group of YSOs can be seen in the IRAC channels, 
but most of them are too deeply embedded to be detected in the 
JHKs image. In fact, only two are visible as very red sources. 
Weaker, but still traceable radio emission is also observed at 
DBlO's position, resulting in an elongation of the contours. 

The paper structure is as follows: in Sect. |2]we describe re- 
duction and photometry for the new JHK$ Sofl images, Sect. [3] 
deals with YSO identification using IRAC photometry, in Sect.|4] 
we describe our Sofl K band spectroscopic observations, which 
are used in Sect.|5]to derive the spectro-photometric distance of 



M. Gennaro et al.: Multiple episodes of star formation in the CN15/16/17 molecular complex. 




Fig. 2. Result of the image filtering. Upper left: K$ band image. 
Upper right: the same image after one allstar pass, with the 
detected sources subtracted. Lower left: The same image after 
spatial filtering. Lower right: the nebulosity that has been filtered 
out in the process. 

CN15/16/17, in Sect. |6] we use sub-mm data from the APE^Q 
Telescope Large Area Survey of the Galaxy (ATLASGAL, 
Schuller et al. 2009) to estimate extinction towards the complex, 
in Sect. |7] we use integrated radio emission from NVSS to put 
additional constraints on the spectral type of the most massive 
members in the DB10 and DB11 clusters, the same is done in 
Sect, fusing the diffuse Bry emission from our spectra, this col- 
lection of data is used to derive the masses of the DB10 and 
DB 1 1 clusters in Sect|9] we discuss the results and summarize 
our conclusions in Sect.fTOl 



2. JHK S Imaging and Photometry 

Observations were performed on the 27th of June 2010. A sum- 
mary of the integration times can be found in Table Q] The inte- 
gration times (DIT) were chosen to be short enough to avoid sat- 
uration of stars brighter than 9 mag in each of the JHK$ bands. 
At each position, NINT=10 frames are averaged to produce a 
single raw frame; NDIT=5 dithering offset positions are used. 
The total exposure time is DITxNDITxNINT 

Given the SofI field of view (FoV) of 4.5'x4.5' and their pro- 
jected distance on the sky of 1.6', both clusters DB10 andDBll 
are covered at each dither position. We spent an equal amount of 
time on the target as on a nearby field located ~ 6' to the East and 
~ 4' to the North from the science field. This control field was 
observed to build an image of the sky free of the nebulous emis- 
sion, which characterizes the CN15/16/17 region. Image reduc- 
tion was performed using a combination of eclipse routines 
(Devillard 2001) and custom-made IDL routines. 



Photometry was performed on the reduced fr ames using th e 
IRAF implementation of the DAOPHOT package (IStetsonlll987h . 
Due to the bright and highly variable background, especially 
in the K$ band, the detection of faint sources in some parts of 
the image was problematic. In order to tackle the problem, we 
performed spatial fi ltering of the image using the Fast-Fourier- 
Transform (FFT, ICoolev & Tukevlll965l) method in IDL. After 
running allstar on the frames once, the task substar au- 
tomatically removes the detected sources producing an image 
which consists of three components: the light from the missed 
sources S(x,y), the nebular diffuse emission N{x,y) and the 
residual noise from the bright infrared sky background S(x,y). 

We used the following procedure to reduce the N{x, y) term. 
In the previous steps we built a model for the stellar point spread 
function (PSF) using bright, isolated sources in the field. The 
FFT of the PSF has been used as a high-pass filter to remove the 
power associated with small wave numbers (i.e. the large-scale 
variations of the diffuse nebulosity). Spatial filtering consists of 
multiplying the FFT of the image with the FFT of the PSF, and 
performing an inverse FFT on this product. In the final image the 
N{x, y) component is significantly reduced. The high-pass filter- 
ing removes most of the power associated with wave numbers 
k < 1 /A, where A is the typical length scale of the PSF (for ex- 
ample the full width at half-maximum (FWHM) in the case of a 
Gaussian PSF). 

The spatial filtering does not account for the S(x,y) noise 
and the Poisson noise associated with the N(x,y) signal. Both 
of these components are indeed highly variable (i.e. their char- 
acteristic length-scale is / < A). The presence of these two 
components sets the fundamental detection limit for the faint 
sources. After applying spatial filtering a second run of the 
daofind peak-finding routine allows detection of additional 
sources, without the need of artificially lowering the detection 
threshold. Only the peak-finding step is run on the processed im- 
ages. The PSF-fitting task allstar is subsequently performed 
on the original, unsubtracted and unfiltered image, using an input 
list obtained by merging the previous detections and the newly 
found peaks. We iterate this procedure one more time to obtain 
the final list of instrumental stellar magnitudes in each band. 

Figure [2] shows the images before and after the spatial fil- 
tering is applied. The upper-left panel shows the original im- 
age in the K$ band. The upper-right panel is the result after one 
PSF-fitting iteration with allstar. The detected sources are fit- 
ted and subtracted from the original image. The residual light 
is coming from the undetected sources, the nebula and the sky 
background. The lower-left panel displays the subtracted image 
after spatial filtering and the lower-right panel shows the sub- 
tracted nebulosity. Schematically we have that: 



Upper Left 



allstar 



Upper right 



->■ Lower left 



APEX is a collaboration between the Max-Planck-Institut fur 
Radioastronomie, the European Southern Observatory and the Onsala 
Space Observatory. 



-< Peak finding -< 

Note that the image cut levels and color-scaling parameters 
in |2] are chosen to be the same in all panels to demonstrate the 
brightness of the subtracted nebulosity and the dimension of its 
luminosity contrast between different areas. 

As already mentioned, our spatial filtering method is not re- 
moving variations in N{x, y) on spatial scales of the order of A. 
Therefore some of the detected sources in a single band might 
be of non stellar origin (i.e. small-scale, large-gradient variations 
in the diffuse nebulosity) which are removed from the final cat- 
alog by matching the three separate lists of detections for J, H 
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Fig. 3. Photometry of CN15/16/17. Black dots are stars identified in J, H and Ks bands, gray dots (upper-right only) are stars 
identified only in H and K$. YSOs are shown in red (class I), green (class II) and yellow (reddened class II). Among the YSOs, the 
circles are sources with JHKs photometry, the squares are objects with only HK$ detections, the diamonds are objects with only 
Ks and the crosses are objects detected neit her in J, H or A"s. The b lack arrows in the near infrared diagrams are the reddening 
vectors for Ak s - 0.5 mag, according to the lRieke & Lebofskvl d 19851) reddening law (suitable for these wavelengths). Upper left: 
Ks vs. J - Ks diagram. Upper right: Ks vs. H - Ks diagram. Lower left: J - H vs. H - Ks diagram. L ower right: I RAC [3.6H4.51 
vs [5.8]-[8.0] color-color diagram. The cyan box indicates th e locus of clas s II objects according to lAllen et al.l (|2P04). The two 
black lines are parallel to the reddening vector derived by the Mathis (1990) e xtinct ion law (suited for mid infrared wavelengths) 
and enclose the region of reddened class II objects according to Megeath et alj d2004h . 



and Ks- For the positional matching we use the Ks band detec- 
tion list as a master catalogue and assume that any real object 
detected in the H band must have a Ks band counterpart and 
any J band detection must have a counterpart in both H and ^s 
bands. A total of 3332 stars are detected in J, H and Ks bands 
and additional 1823 stars are detected in H and Ks only. 

Instrumental magnitudes have been transformed into the 
2MASS photometric system by matching our detections with 
bright and isolated 2MASS sources. A linear fit of the 
2MASS vs. instrumental magnitudes was performed between the 
matched stars, allowing for a color term. The resulting photom- 
etry for the entire region is displayed in Fig. [3] 



Table 1. Observing log for the JHKs observations. 



Band Exposure Time Seeing 
[s] [arcsec] 



H 



500 
300 

125 
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3. IRAC mid infrared photometry 

In addition to the near infrared photometry from SofI, 
analysed the mid infrared photometry for the region in the 
IRAC bands centered at 3.6 |am, 4.5|im, 5.8|om and 8. Op 
The region has been observed within the GL IMPSE II sun 
dBenjamin etai]|2003l:lc"hurchwell et aTll2009l) . In this work, 
used the GLIMPSE II Epoch 2 November '09 Point Sou 
Catalog. 

In an area of 6' x 6', centered on DB11, 342 sources are i 
tected in all of the 4 IRAC bands. Among these objects, clas 
and class II YSOs can be classified using their mid infrared ( 
cess due to the circumstellar envelopes and disks, respective 
lAllen et all (l2004h showed that these two classes of objects i 
cupy different areas in the [3.6]-[4.5] vs. [5.8]-[8.0] color-co 
diagram with respect to "purely photospheric" objects which 
clustered around the (0,0) position. Even though there mij 
be a small overlap b etween the class I and class II obje< 
Megeathet al. (2004) proposed a working criterion to ident 
the two classes of objects and separate them from the mi 
evolved, envelope and disk-free stellar sourc es. Using the m( 
els of the emission of class I and II objects bv l Allen et alj d200 I 
they identify a rectangular area with 0.0 mag < [3.6] - [4.5^ 
0.8 mag and 0.4 mag < [5.8] - [8.0] < 1.1 mag as the locus 
class II YSOs. This box is displayed in the lower-right panel ui 
Fig. [3] The area of the diagram above and to the right of the box 
is occupied by class I objects. Part of the class II objects resides 
in a reddening stripe on the upper left of the unreddened class II 
box. The two parallel lines in t he same figure show the direction 
of the reddening vector for the Mathis (1990) extinction law. 

We identify 13 class I, 9 class II and 6 reddened class II ob- 
jects in our sample. Class I: 2 sources have JHKs photometry, 4 
HK$ photometry and a further two are only detected in Ks band. 
Class II: of the 9 class II sources 7 have JHKs photometry and 2 
have HKs photometry. Reddened class II: these object have very 
large visual extinction values. A reddening of 0.4 mag in [3.6]- 
[4.5] corresponds to about 30 mag in Ay. The average color of 
these reddened sources is [3.6]-[4.5] ~ 1.1 mag, while for the 
unreddened class II objects an average [3.6]-[4.5] ~ 0.4 mag 
is observed. This typical ~ 0.7 mag reddening corresponds to 
Aj sa 14.8 mag, Ah sa 9 .2 mag and Ak s ~ 5.9 mag, using the 
Riek e & Lebof skv (1985) extinction law. Such large extinction 
values hamper the detection of reddened class II objects in our 
SofI observations: one of them is detected in JHKs, while 3 are 
detected in Ks band only. 

Given that the region is projected close to the Galactic Center 
(/ = 0?58, b - -0?85), some contamination from asymptotic gi- 
ant branch (AGB) star s might be present in th e photometric YSO 
sample. According to Maren go et al.l (|2008), the loci of oxygen 
rich AGB stars and supergiants slightly overlap with class II 
YSOs. Only spectroscopy of these sources could help discrim- 
inate between different stellar types. The class I sample, on the 
other hand, should not be contaminated. 

The IRAC 8(a.m image of the region is presented in Fig. |4] 
with the positions of the identified YSOs superimposed. Their 
spatial distribution shows that star formation is occurring 
throughout the entire region and not only in the DB 1 and DB 1 1 
clusters. For example, a small clustering of YSOs is visible to- 
wards the southern part of the region. This area appears very dark 
in the near infrared images, due to extinction from the parental 



2 Available at the following URLs: 
irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/ (images) 
irsa.ipac.caltech.edu/cgi-bin/Gator/nph- 
scan?submit=Select&projshort=SPITZER (catalogs) 
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Fig. 4. IRAC 8 \irs\ image with the marked positions of the iden- 
tified YSOs. Crosses correspond to objects which are not iden- 
tified in the JHKs bands. Diamonds are objects identified only 
in Ks, squares are objects identified in H and Ks, circles are ob- 
jects identified in J, H and Ks- Class I objects are in red, class 
II in green, reddened class II sources in yellow. The dashed box 
corresponds to the observed SofI field. The contours trace the 
emission of cold dust at 870 [im from the ATLASGAL survey 
of the Galactic Plane. Contour levels are between 0.5 and 9.5 
Jy/beam in steps of 1 .0 Jy/beam. 



molecular cloud (see Sect. |6}. The DB11 and DB10 clusters 
show a paucity of YSOs. Part of the reason for this could be 
that the two clusters are evolved and effectively not hosting any 
YSO. However, the lack of YSOs might be partially due to the 
high brightness of the nebula in the mid infrared at the positions 
of DB10 and DB11, which hampers source detection. 



4. K band spectroscopy 

Spectroscopic observations with SofI were performed on the 
27th of June 2010. The typical seeing during the observations 
varied between l'/2 and l'/6. The spectroscopic targets were pre- 
viously selected from the 2MASS point source catalog by choos- 
ing the brightest stars at the DB10 and DB11 cluster locations. 
We aimed to observe massive stars in the CN15/16/17 region 
and, by classifying their spectral type, infer the region's distance 
as well as the extinction towards those objects. To minimize the 
contamination by foreground main sequence stars, only candi- 
dates with H - Ks > 0.3 mag were considered. However, it was 
still possible for field giants to contaminate the sample, as they 
have similar photometric properties as more distant and more 
reddened massive main sequence stars. Five different slit posi- 
tions were used; the positions were chosen in order to observe 
the 10 brightest stars in the region, but given the high degree 
of crowding more stars fell into the slits and we could extract 
spectra for 3 additional stars with sufficient signal-to-noise ratio 
(SNR) to allow for spectral classification. 
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RA (J2000) 

Fig. 5. K& band image of CN15/16/17. Markers show the posi- 
tions of the stars with spectral classifications. Green crosses are 
early type stars, yellow circles are candidate YSOs, red squares 
are red giants, and the magenta diamond is an AGB star. IDs 
correspond to the first column of Table [2] The two white dashed 
circular regions are the areas corresponding to clusters DB10 
(lower right, 21" in radius, i.e. 0.12 pc at 1200 pc distance) and 
DB1 1 (center, 32" in radius, i.e. 0.18 pc at 1200 pc distance). 



We used the spectroscopic mode of NTT/SofI selecting the 
third order of the HR Grism, with a dispersion of 4.63 A/pixel in 
the 2.00-2. 30|im wavelength range. We always used an 076 slit, 
resulting in a spectral resolution of 2200. Several lines that allow 
for spectral classification of early type stars were present in the 
chosen spectral range, such as Bry (H i, 2.166|im), and He i 
lines (2.059, 2.1 13(xm). Other lines which were present allowed 
for the classification of the late type giants contaminants, such 
as Na i (2.206, 2.209 |xm) and Ca i (2.261, 2.263, 2.265 urn), as 
well as the CO absorption band at 2.3ii.m. 

For each slit multiple nodding cycles (NC) were observed. 
The slit was first put in position A and then shifted along its 
long axis to position B. Removal of the sky background was 
performed using the difference between the A and B positions. 
The number of NC was optimized, depending on the brightness 
of the sources, to obtain a SNR sufficient to perform spectral 
classification. 

Reduction wa s performed using the ESO gasgano pipeline 
dlzzo et al.l l2004h . The pipeline applied flat fielding and dark 
subtraction for each raw image, and used xenon and neon arc 
lamp images to correct for slit curvature. For each nodding cy- 
cle the (A-B) + (B-A) sky-subtracted images were obtained. The 
final step was to shift and sum of the 2 x NC sky-subtracted im- 
ages. Wavelength calibration was also done using the xenon and 
neon arc lamp spectra. To extract the spectra we used the IDL - 
optspecextr optimal extraction packagqfj. This software traces 



3 The software is publicly available at: 
http://physics.ucf.edU/~jh/ast/software/optspecextr-0.3.l/doc/ 



the sp ectra and extracts the signal using the algorithm of lHornel 
( 1986) to optimize the SNR ratio of the final spectrum. 

Throughout the observing run we obtained spectra of late B 
and early A type reference stars in order to correct for telluric ab- 
sorption. These stars have the advantage of an almost featureless 
spectrum in the K band. The positions of the standard stars were 
chosen to minimize the airmass difference with the CN15/16/17 
targets. Telluric standards were observed about every 30 min- 
utes in order to have similar atmospheric conditions between the 
standards and the targets. The only prominent feature in these 
early type telluric standards is their Bry photospheric absorp- 
tion. Therefore the extracted standard star spectra were first cor- 
rected by removing the Bry line using the IRAF/splot task to 
fit a Voigt profile to the Bry line and subtract the result of the 
fit from the spectrum. We then used the IRAF/telluric task to 
optimize the telluric correction. The spectrum of each target star 
was divided by the telluric standard star spectrum (with Bry sub- 
tracted). The telluric task allows iterative optimization of the 
atmospheric correction to compensate for slight variations in air- 
mass and atmospheric conditions between the science target and 
the standard star observation. This is accomplished by scaling 
the amount of atmospheric absorption and shifting the spectrum 
in wavelength until a satisfactory correction is achieved. 

The SNR per spectral resolution element required to achieve 
a good spectral classification was of the order of 100. Due to 
bad seeing conditions throughout the observing run, the achieved 
SNR is lower, ranging between 50 and 80 (depending on the 
stars). This led to a typical uncertainty of ~ 2 spectral subclasses 
for the classification of the early type stars in our sample. 

5. Spectral Classification: distance and extinction 
to CN1 5/1 6/1 7 

The stars for which spectral classification was possible are 
marked in Fig. [5] superimposed on the SofI Ks band image of the 
region. The extracted and corrected spectra are shown in Fig. [6] 
A total of four early B-type stars were identified and four addi- 
tional YSOs candidates. The rest of the classified stars are back- 
ground giants, four classified as red giants and one as an AGB 
star. The results of the classification procedure are summarized 
in Table [2] 

Classification of the early type stars was done by comparing 
the observed spectra with the medium resolution near infrared 
spectral atlas of O and early B Stars by Hans on et al.l (|2005), 
complemented by spectra from Bik et al. (2005J). To compare 
these high quality spectra with our lower SNR, lower resolution 
ones we used the following procedure: we first convolved the at- 
las spectra with a Gaussian kernel with a <x equal to 4.6 A/pixel, 
i.e. the SofI wavelength dispersion in K band. Second we mea- 
sured the SNR in a featureless part of the observed spectrum and 
added noise to the atlas spectra until they reached a comparable 
SNR. 

The spectral classification was then performed by visual 
comparison of the observed spectra with these degraded atlas 
spectra. The strength of the Bry absorption line and of sev- 
eral helium lines in the observed spectral range (He i 2.059ixm, 
2.113|im) was used to get the best match. The Bry absorption 
lines of the early type stars inside DB10 and DB1 1 are contam- 
inated by the nebular Bry emission from hydrogen recombina- 
tion. This emission, though, is quite narrow and does not affect 
the wings of the photospheric absorption Bry lines. Hence the 
absorption win gs can still be used for spectral classification (see 
lBiketal.ll2012l for a similar classification procedure). An ex- 
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Fig. 6. Spectra of the classified objects. The IDs correspond to 
the first column of Table [2] 
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Fig. 7. The Bry line spectral region for the El source compared 
to the degraded spectrum of a BOV star template. 



iLejeune & Schaererl (1200 ll) ZAMS of solar metallicity. The er- 
ror in the estimated extinction and distance mainly comes from 
the uncertainty in the derived spectral type. 

We estimate an average distance to the region using the 
spectro-photometric values from our spectral classification. We 
used the reasonable assumption that the four stars are all located 
at the same distance from the Sun. By considering the central 
value of the distance interval for each of the four classified early 
B stars we obtained d — 1.2 + 0.5kpc. The quoted uncertainty is 
the standard deviation of the four central values. The total visual 
extinction for the only B star identified in DB 1 1 (El)is Ay - 4.2 
mag, while the three B stars identified in DB 10 (E2, E3, E4) pro- 
vide a slightly larger average value of the extinction of Ay ~ 5.5 
mag, consistent with DB10 being younger and more embedded. 

The identified YSOs are characterized by featureless, red 
spectra (see Fig. [6}. Y2, Y3, Y4 show Bry emission and Y2 also 
shows He I emission at 2.059 and 2.1 13|J.m. This emission might 
come from the nebula, but might also come from the circumstel- 
lar material, with a similar mechanism as observed in Be stars 
(Clark & Steele 2000). 

For the classifi cation of giants we used the atlas by 
Ravn er et al.l (|2Q09). Given the similar resolution between our 
observations (R = 2200) and the spectra of this atlas (R = 2000), 
we could classify them by direct comparison of the observable 
features, such as the CO absorption band at 2.3 |xm, the Na i 
(2.206, 2.209 urn) and the Ca i (2.261, 2.263, 2.265 urn) absorp- 
tion lines. 



ample of such comparison of the Bry absorption line wings is 
shown in Fig. [7] where we show the Bry line for the El source, 
classified as a B0V-B2V star, and the degraded atlas spectrum 
for a BOV type star. 

Given the spectral types derived for the early type stars, it 
is possible to obtain an estimate of their distances and individual 
reddening values. Using the intrinsic colors of early type stars by 
iKenvon & Hartmanri (Q995) we d etermined the color exces s of 
the observed stars. Adopting the Rieke & Lebofsky ( 1985J) ex- 
tinction law we then estimate the foreground extinction. Given 
the low values of the color excess (see Table |2j, the choice of 
the reddening law was not critical, and the differences expected 
from different extinction laws are small compared to the other 
sources of uncertainty. The distance was derived by the differ- 
ence between the dereddened magnitude and the absolute mag- 
nitude expected for the given spectral type. For the latter we used 
Leieune & Schaerer (2001) models. We checked that the results 
did not depend on the actual value of the age of the isochrone 
adopted, in the range 1 to 3 Myr. Results in Table [2] are for the 



6. Sub-mm continuum emission 

The cold dust emissio n at 870 yim from the ATLASGAL survey 
of the Galactic Plane ( Schull er et al.l|20 09) can be used to trace 
extinction. The 870(a.m dust emission from the region is over- 
plotted in Fig. [4] The cold dust traces the positions of the identi- 
fied YSOs in the southern part of the CN15/16/17 complex very 
well. Also the more isolated YSOs in the northern area are sur- 
rounded by emission at 870 |xm. The peak of 870 \xm emission 
is coincident with the bright Ultra Compact H n region in the 
south-west of the field. This dense material is therefore associ- 
ated to a region of high mass star formation. The remaining YSO 
population in the south of the field follows the arc-like southern 
region of 870 |xm emission and is not associated with H n emis- 
sion. This is visible in Fig. Q] where no radio emission can be 
observed in this area. Hence these YSOs are likely less massive 
objects, and are not able to ionize their surroundings. 

Using equation (1) of Schuller et al. (2009) with the de- 
rived spectro-photometric distance of 1.2 kpc and assuming 
k v - 1.85cm 2 g _1 , we estimate a mass of w 1O5M0 for the south- 
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Table 2. Spectral Types of the classified stars. For early type stars the inferred extinction and distances are indicated as well. 



ID 



R. A. (2000) 

h m s 



Dec.(2000) 



Spectral Type 



[mag] 



J-K s 
[mag] 



H-K s 
[mag] 



Ay 

[mag] 



i£ s 



[mag] 



Distance 
[kpc] 



Early Type 



El 17:50:24.03 


-28:53:07.8 




B0V - B2V 




10.36 


0.73 


0.21 


4.2 ±0.1 


E2 17:50:17.53 


-28:53:39.2 




B1V-B3V 




10.44 


0.87 


0.32 


5.7 + 0.1 


E3 17:50:16.94 


-28:53:40.4 




B1V-B3V 




9.75 


0.93 


0.26 


4.7 + 0.1 


E4 17:50:16.97 


-28:53:41.9 




B5V-B9V 




11.08 


0.89 


0.36 


6.0 + 0.1 


Young Stellar Objects 


















Yl 17:50:18.68 


-28:53:31.8 




red spectrum 




10.82 


4.63 


1.51 


- 


Y2 17:50:24.35 


-28:53:03.1 


Bry, 


He i 2.059, 2.113 
(all in emission) 


urn 


11.70 


2.46 


1.31 


- 


Y3 17:50:23.94 


-28:53:06.1 


featureless, Bry emission 


13.49 


1.67 


0.96 


- 


Y4 17:50:15.41 


-28:54:11.9 


featureless, Bry emission 


11.48 


0.97 


0.36 


- 


Asymptotic giant branch 


. 
















AGB1 17:50:13.75 


-28:51:38.0 




M5 - M9 III 




8.86 


5.35 


1.80 


- 


Red Giant Branch 


















RG1 17:50:21.94 


-28:53:06.3 




K2 - K4 III 




9.73 


4.56 


1.49 


_ 


RG2 17:50:21.90 


-28:52:58.5 




K2 - K4 III 




9.99 


4.82 


1.78 


- 


RG3 17:50:17.17 


-28:52:37.0 




K2 - K4 III 




9.98 


4.73 


1.55 


- 


RG4 17:50:22.48 


-28:53:18.0 




K5 - K8 III 




11.55 


3.53 


1.18 


- 



0.45 ±0.01 1.3-2.4 

0.62 ±0.02 1.0-1.6 

0.51 ±0.01 0.7-1.2 

0.65 ±0.02 0.6-1.3 



west clump associated with the Ultra Compact H n region. The 
mass of the whole arc is estimated to be as 230M© by integrating 
the flux over the region for which the intensity is larger than 2 
Jy/beam. Given the presence of nearby hot stars, we assumed a 
dust temperature of 50K, keeping all the other parameters as in 
Schu ller et alj d2009l) . Such a large temperature for the cold gas 
comp onent is typical of r egions of massive star formation (see 
e.g. lSridharan et al.ll2002l) . 



Using equation (2) of ISchuller et al.1 (120091) . we additionally 
derived the H2 column density towards the region, and then used 
N(U2)/A V = 0.94 x lO^cm^mag- 1 (Frerkin getai]|1982l) to 
convert it into visual extinction. At the position of the El star 
we obtain Ay ~ 4.2 mag, in perfect agreement with the spec- 
troscopic measurements. In comparison, the ATLASGAL maps 
provide very high values for the extinction along the lines of 
sight towards E2, E3 and E4 in the DB10 cluster region, with 
Ay ~ 15 mag. The extinction values depend on the choice of 
the dust temperature. For example, choosing T - 30 K, the esti- 
mated Ays are almost doubled. Therefore our quoted extinction 
values from sub-mm emission have to be considered as approx- 
imate estimates. The fact that in DB11 the spectroscopic and 
sub-mm results are in agreement, while they are not for DB10, 
might be due to several factors. For example the bright nebula 
might partially affect our photometric measurements and there- 
fore the reddening estimates. However, a difference of 10 mag in 
Ay would correspond to an error in the J-Ks color of more than 
2 mag, which is too large for these bright sources. Another pos- 
sibility is that the differences in the extinction derived for DB10 
are related to the beam size of the ATLASGAL maps of 19". 
This is comparable in size to the visual radius of the DB10 clus- 
ter. Given this comparably large beam, the measurements of the 
flux at the DB10 position are possibly affected by the emission 
of the prominent nearby south-west arc and clump. DB11, on 
the contrary is located in a region of lower sub-mm brightness, 
with no close prominent features, hence the spectro-photometric 
results and those from the 870 \xm maps are very similar. 



7. Integrated radio continuum flux 

The radio flux from the two H 11 regions in the CN15/16/17 com- 
plex can be used to put additional constraints on the spectral type 
of the most luminous ionizing source. The amou nt of ionizing 
flux is a very steep function of the spectral type dMartins et al.l 
2005). Consequently, the total ionizing flux within an H 11 region 
is dominated by the hottest, earliest type star. 

The free electrons inside the H 11 region interact with the ions 
producing Bremsstrahlung. This free-free radiation flux, mea- 
sured at a given wavelength, can be converted into an ionizing 
flux. Assuming that each energetic photon ionizes an atom (no 
leakage), that the nebula is optically thin for free-free contin- 
uum emission and that ionization and recombination are in equi- 
librium, the amount of energetic photons and radio continuum 
emission can be related to each other. 

The values of the in tegrated radio flux a t 1 .4 GHz are taken 
from the NVSS catalog (ICondonet al.lll998l) . The ultra compact 
H 11 region is the brightest at 1 .4 GHz, with a flux of 657 mJy, 
while the flux of the H 11 region associated with DB11 is 572 
mJy. The observed radio flux has first to be converted into an 
emitted flux by multiplying it with And 2 . 

We used the relations bv lKurtz etal.1 (Q~994) to obtain an es- 
timate of the ionizing fluxes. Given the estimated distance of 
d — 1.2 + 0.5kpc, (see Sect. [5]) we obtain: 

log Go BU [photons s" 1 ] = 46.8^J 



+0.3 



log Go ""[photons s" 1 ] = 46.91^ 

where the flux Qo is the number of photons with A < 912 A emit- 
ted per second. 

We compared these values with those predicted by theo- 
retical models of hot stars and derived estimates for the spec- 
tral types of the ionizing sources inside the DB 1 1 cluster and 
the Ultra Compact H 11 region. Adopting the calculations by 
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iMartins et al.l (120051) . the expected flux from the latest spec- 
tral type available -an 09.5V- varies between log Qq ~ 47.56 
and \ogQo ~ 47.88, depending on the adopted T e g scale (de- 
rived either from models or observati ons). Similarly we con- 
sidere d the ionizing flux computed by Diaz-Mil ler et al.l (1998) 
using Kurucz (1993) model atmospheres. These computations 
extend to lower stellar temperatures than those of Martins et al. 
(2005). The authors predict an ionizing flux log Qq — 47.02 for 
a 17.5Mq, solar metallicity, main sequence star, corresponding 
to r e ff = 30, 000 K. This flux is consistent with our estimates. 
No spectral type is quoted by the authors, but from [Martins et al] 
d2005l) a r eff = 30, 000 K corresponds to an 09.5 class V star, for 
which, however, their predicted inonizing flux is log go ~ 47.56. 
These differences testify the difficulties in modeling the atmo- 
spheres of hot stars. Given this uncertainty we cannot constrain 
the expected spectral type to better than few spectral subclasses, 
but we can safely assume that our estimated values are consis- 
tent with the earliest type of the classified stars being late O or 
early B. 



8. Nebular emission 

The intensity of the diffuse nebular emission can be used to 
further constrain the spectral t ype for the br i ghtest stars inside 
DB 10 and DB 1 1 . As shown bv lHanson eTall 00021) . the ratio of 
the fluxes of the Bry and He i 2.1 13ii.m lines is an indicator of 
the temperature of the ionizing source. 



We measured the line fluxes in the nebulae of DB 1 1 and 
DB10 using the parts of the 2D long-slit spectra which spatially 
coincide with the nebular regions. In order to obtain the maxi- 
mum SNR and thus the best possible constraints on these two 
fluxes, we combined all five slit positions for which we obtained 
spectra. Since exposure times and observing conditions were dif- 
ferent among the five slit positions, we first calibrated the 2D 
frames in flux. We used the spectra of the telluric standard that 
was observed immediately after or before the given slit. We inte- 
grated the measured counts convolved with the 2MASS Ks filter 
response curve and used the available magnitude for the star to 
obtain an absolute calibration of the flux for the standard star. By 
dividing the 2D frames for this calibrated standard spectrum and 
adjusting for the different exposure times we obtained 2D flux 
calibrated frames. We identified in each 2D frame the regions of 
diffuse Bry emission associated with either DB10 or DB1 1, and 
summed the flux in each region. The fluxes for the five slits were 
then co-added to obtain the total nebular flux. The summed spec- 
tra forDBlO and DB11 are shown in the top panels of Fig. [8] 

According to lHanson et~ai1 d2002l) a He i (2. 1 1 3)/Bry flux ra- 
tio between 1% and 3% implies that some late O type stars con- 
tribute to the ionization of H n regions, and lower ratios imply 
that no O star is present. In both DB 1 and DB 1 1 , we find no He i 
(2.1 13) emission, shown in the lower panels of Fig. [8] Assuming 
a similar line shape and FWHM to the Br-gamma line, we plot 
the expected flux level for the 1% and 3% flux ratio cases. We 
also show the lcr and 3<x detection levels, where the standard 
deviation <x of the flux is measured in adjacent featureless parts 
of the spectrum. In both DB10 and DB11, a 3% flux ratio can 
be excluded at the l-<x level, but a 1% flux ratio cannot be ruled 
out by our data. We can therefore exclude the presence of early 
O type stars as the ionizing source in these regions. The earliest 
spectral type would be early B or late O-type stars. 



9. Mass estimates of the DB10 and DB11 clusters 

In this section we will focus on the two clusters in the region, 
DB10 -more embedded and compact- and DB11 -less embed- 
ded, more extended. We consider two circular regions as the 
clusters areas. These areas correspond to the diffuse nebulosity 
areas visible in the Ks band image (see Fig.|5j- DB10 is centered 
at (RA, Dec)j2ooo = (17 h 50 m 17. s 64, -28°53'41'/6), (/,fc)gai = 
(0.564,-0.854), with a visual radius of 21". Object DB11 is 
centered at (RA, Dec) J20 oo = (17 h 50 m 23. s 79, -28°53'05'.'l), 
(/,fc)gai = (0.584, -0.868), with a visual radius of 32". At a dis- 
tance of 1200 pc the radii of the clusters are ~ 0.12 pc and ~ 0.18 
pc for DB10 and DB1 1 respectively. 

The near infrared color-magnitude and color-color diagrams 
for the two cl usters are shown in Fig. |9] Overplotted is a 1 Myr 
isochrone bv iLejeune & Schaererl (12001 1) . We have chosen this 
age based on the presence of the molecular material in the sur- 
rounding bubble. This material will be removed by the winds 
and radiation pressure from massive stars. A slightly older age, 
up to few Myr might still be possible but we consider 1 Myr as a 
typical value. A more precise age determination can be obtained 
with deeper and higher angular resolution images revealing the 
pre main sequence population. We adopt the average value of 
the spectro-photometric distance obtained from the classified B 
stars, d - 1 .2 kpc. Using the same the spectro-photometric anal- 
ysis and the Rieke & Lebofskv ( 1985) extinction law we derive 
the extinction towards the clusters. For DB10 we used the aver- 
age extinction of the stars E2, E3, E4, Ay - 5.5 mag (Aj - 1.54 
mag, Ah - 0.96 mag and A# s = 0.61 mag). For DB11 we 
adopted the value from El, Ay - 4.2 mag (Ay = 1.18 mag, 
A H - 0.74 mag and A Ks - 0.47 mag). 

The positions of the black squares in Fig.|9]indicate the mag- 
nitudes of the classified objects, with early type stars in black, 
YSOs in green and red giants in red. The bright nebulosity in 
DB10 hampers the detection of faint Ks band sources, conse- 
quently the diagrams appear quite empty in redder areas. In 
DB11 the nebulosity is not as bright and therefore fainter (and 
redder) sources can be seen in the diagram. 

The classified B stars nicely follow the 1 Myr isochrone. The 
spectroscopically identified YSOs appear redder than the main 
sequence locus in both the Ks vs J - Ks and the Ks vs H - Ks di- 
agrams. From the color-color diagram, in the case of DB 1 1 , it is 
clear that the redder YSO colors are not due to extinction but are 
related to near infrared excess from the circumstellar material. In 
this diagram the two DB 1 1 YSOs (Y2, Y3) are indeed located 
below the reddening sequence. The YSO Y2 is identified in the 
GLIMPSE II catalog as well. With its colors of [5.8]-[8.0] = 1.01 
mag and [3.6]-[4.5] = 0.69 mag it is one of the objects identified 
as a class II YSO in the IRAC color-color diagram of Fig. [3] 
In the DB1 1 diagrams of Fig.|9]Y2 is the redder (J - K s = 2.47 
mag) of the two marked YSOs. The YSO in DB10 (Yl) does not 
show the same near infrared excess. This object might indeed 
be more embedded as it appears from the high extinction in the 
color-color diagram. Unfortunately, we could not find a match 
for this source in the GLIMPSE II catalog, but the very red and 
featureless spectrum indicates a class I type. The three red gi- 
ants classified in the DB1 1 area (RG1, RG2, RG4) all appear to 
be very reddened background sources, most likely belonging to 
the Galactic Bulge. 

The cluster DB1 1 was already studied in lDutra et al.l (12003) 
without the help of spectroscopy and with very different results. 
After performing statistical field subtraction, the authors iden- 
tified a much redder main sequence for the region with a color 
of H - Ks ~ 1.12 mag. Fitting this sequence they derived an 
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Fig. 8. Spectra of the diffuse nebulae in the DB10 (left) and DB1 1 (right) regions. The positions of the Bry and the He i 2.1 13 lines 
are marked by vertical dotted lines. The lower panels are a zoomed-in version in the wavelength region of the He i 2. 1 13 line. The 
horizontal dashed lines represent 1 and 3<x levels above the continuum. The squares and the circles represent the expected peak 
values of the He i 2.113 line in the case of a 1% and 3% He i/Bry flux ratio, respectively. 



Ay of 15 mag (A^ s = 1.68 mag). The classified B star in DB11 
(El) shows, in contrast, a much bluer H - Ks color of only 0.21 
mag and a much lower extinction Ay - 4.6 mag (A^ s = 0.47). 
Three of the classified objects lie in the region H - Ks ~ 1.12. 
One of them is a red giant and the other two are the Y2 and 
Y3 YSOs. The presence of the first is an indication that this part 
of the DB 1 1 color magnitude diagram is contaminated by red gi- 
ants. Moreover the color-color diagram shows that the two YSOs 
clearly have a Ks band excess, being displaced from the main se- 
quence by ~ 0.6 mag. If the main sequence of DB1 1 in the Ks 
vs H -Ks diagram were indeed at H - Ks ~ 1.12 mag we would 
expect the two YSOs to be ~ 0.6 mag redder than they are. 

Assuming that the 10th brightest star of the group that they 
identified as cluster members was a B0V star. lDutra et al.l (12003) 
derived a distance to the cluster of 7.6 kpc. This scenario was 
corroborated by the valu e of the kin e matic distance for the asso- 
ciated H ii region. iKuchar & Clarkl (1 19971) estimated a velocity 
along the line of sight of 15 km/s. This value of the velocity 
can be translated into a distance once a model of the Galactic 
Rotation curve is adopted. Given the line of sight for this region 
-very close to the Galactic Center- a very small radial velocity is 
expected for all orbits but those very close to the Galactic center. 
The derived kinematic distances are indeed 7.5 kpc (near) and 
9.5 kp c (far) using the Galactic rotation model by iBrand & Blitz! 
(Il993l) with a distance of the Sun from the Galactic Center 
Ro = 8.5 kpc and an orbital velocity at Ro, 8 = 220 km s -1 . 
The kinematic distance value is very uncertain though. For ob- 
jects that are moving almost perpendicular to the line of sight, 
the peculiar orbital motions and the internal motion of the H n 
region can dominate the value of the line of sight velocity. 



Our spectroscopic observations give contrasting results for 
the distan ce of the clu s ter. W e argue that the observed main se- 
quence of lDutra et al.l (|2003) might be a result of the subtrac- 
tion technique, related to the difficulties of choosing a represen- 
tative control field, when looking towards the Galactic Center. 
Along such lines of sight the extinction pattern might be very 
patchy and strong differences can be observed between adjacent 
areas. This is also the case for the CN15/16/17 region where, as 
visible in Fig. [4] local concentration of cold dust correspond to 
large values of extinction. The area centered on DB1 1 has been 
cleared of the molecular material by the early type stars with 
their winds. This material accumulates in the rims of the bubble, 
but not directly at the cluster position; as a consequence extinc- 
tion is smaller in the cluster's area. Therefore it might be pos- 
sible that a larger fraction of background stars is visible along 
the DB11 line of sight, when compared to the neighboring re- 
gions used by Dutra et al. (2003) for statistical field subtraction 
(see also Fig. Q}. This would result in an under-subtraction of 
background giants in the DB1 1 color-magnitude diagram. 

If the 10-th most massive s tar of the cluster w ould have been 
of B0V type -as proposed bv lDutraet al. (2003)- the presence 
of several O stars would be expected. We ruled this out by show- 
ing that both the nebular emission and the radio flux are consis- 
tent with the brightest stars being B-type. In our target selection 
for spectroscopy we have chosen the brightest stars in the re- 
gion and identified them as early B stars. Nevertheless we might 
have missed some earlier-type stars, because they could be more 
extincted and appear fainter. However, in this case the diffuse 
emission and the radio data would not have been consistent with 
DB1 1 hosting at most early B or very late O stars and being at a 
distance of 1 .2 kpc. A distance of 7.5 kpc would imply the pres- 
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ence of early O stars, which would have been detected through a 
much larger He i nebular emission at 2.113|^m (see Sect. [8}. 

Assuming that we identified the most massive star in DB10 
(B1V, ~ 12M©) and DB11 (BOV, ~ 15M ) and using the rela- 
tion between the mass of the most massive star in a cluster and its 
total mass bv iLarsonl (120031) . we can give an estimate of the clus- 
ters masses. We obtain M DB io ~ 170M© and M DB n ~ 275M©, 
both comparable to the mass of the <x Ori cluster in the Orion 
OB lb association, Mo-ori ~ 225 + 30M Q (Sherr y et al]|2004l) . 



10. Discussion and Conclusions 

We performed a multiple wavelength study of the young star 
forming region CN15/16/17 towards the Galactic Center. We de- 
rived a consistent picture of the region combining near infrared 



spectroscopy, near and mid infrared photometry, sub-mm con- 
tinuum emission and integrated radio fluxes. 

The region hosts two near infrared clusters (DB10 and 
DB11, visible in the Sofl images) and a loose association of 
YSOs identified in the GLIMPSE II point-sources catalog. The 
DB11 cluster is associated with an H n region visible in the 
NVSS radio continuum maps. Another H n region is also vis- 
ible, associated to a small group of very embedded mid infrared 
bright sources. From ATLASGAL sub-mm contours we identify 
an arc-like structure in the south-western part of the region. This 
arc of cold, dense material is associated with the southern group 
of YSOs. 

Considering this collection of evidence from multiple wave- 
length regimes, a qualitative picture of the star formation history 
of the region can be traced. The central cluster, DB 1 1 , is the old- 
est star formation site, followed by the younger cluster DB10 
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and the youngest Ultra Compact H n region. Star formation is 
proceeding sequentially from DB11 outwards. This is corrobo- 
rated by the presence of the cold dust arc in the south, which 
is associated, in turn, with the most populous group of YSOs 
of the region. In addition, two small sub-groups of YSOs are 
visible in the north and north-west, again associated with cold 
dust emission. These two smaller concentrations are at a similar 
projected distance from DB11 as the south-west arc. However, 
the age determination of the clusters is quite uncertain and, e.g. 
spectroscopic identification of pre-main sequence members is 
required to place DB11 and DB10 on a more robust age scale. 
Further investigation is required to understand whether the ap- 
parently | _y_ojHig^r_eTCnts_Jwv£_be£njrt£g^r^ by_fhe older ones 
(see e.g lWatson et alj|201ft iDeharven g et al.ll2010l as examples 
of analysis of triggered star formation in mid infrared bubbles), 
or star formation has been seeded independently in multiple sites 
of the parental molecular cloud, in analogy, e.g., to w hat is ob- 
served in the Rosette Mo lecular Complex (Li & Smith 2005; 
iRoman-Zuniga et al.ll2008h . 

Through spectroscopic characterization of the brightest 
members of the DB10 and DB1 1 clusters, we obtained a spectro- 
photometric distance to the region as well as the extinction value 
towards the clusters. We found that the region is at d — 1.2 + 0.5 
kpc, i.e. much closer than previously thought when only H and 
Ks photometry was used for the analysis of the DB 1 1 cluster. 
From the mass of the most massive stars we estimate the stellar 
cluster masses to be Mdbio ~ 170Mq and Mdbh ~ 275Mq. 

This work affirms the importance of the multi-wavelength 
approach to studies of young star forming regions. We illustrate 
that the combination of imaging and spectroscopy is essential to 
the understanding of embedded stellar populations. 
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